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An Alternative to the e n Method for Determining 

Onset of Transition 

E. S. Warren* and H. A. Hassan t 
North Carolina State University, Raleigh, NC 27695-7910 


Abstract 

A new approach has been developed to determine 
the onset of transition. The approach is based on 
a two-equation model similar to those used in the 
study of turbulence. The approach incorporates in- 
formation from linear stability theory on streamwise 
or Tollmien-Schlichting instabilities. The present 
approach has proven to be an inexpensive alterna- 
tive to the e n method for determining the onset 
of transition on a flat plate and airfoil for a vari- 
ety of Reynolds numbers and freest ream intensities. 
Further, the method is incorporated into two flow 
solvers, boundary layer and Navier-Stokes. This 
made it possible to calculate the laminar, transi- 
tional, and turbulent regions in a single computa- 
tion. 

Introduction 

The process of transition from laminar to turbu- 
lent flow remains one of the most important unsolved 
problems in fluid mechanics and aerodynamics. Vir- 
tually ail flows of engineering interest transition from 
laminar to turbulent flow. Because transitioned flows 
are characterized by increased skin friction and heat 
transfer, the accurate determination of heating rates 
and drag critically depend upon the ability to pre- 
dict the onset and extent of transition. However, 
no mathematical model exists which can accurately 
predict the location of transition under a wide range 
of conditions. Design engineers resort to methods 
which are based on either empirical correlations or 
linear stability theory. 

The e n method is currently the method of choice 
for determining transition onset. The method 
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is based on linear stability theory and generally 
requires the following steps 

(i) Pre-calculation of mean flow at a large number 
of streamwise locations along the body of inter- 
est. 

(ii) At each streamwise station, a local linear sta- 
bility analysis is performed. By assumptions 
of the linear theory, the unsteady disturbances 
are decomposed into separate normal modes of 
different frequency. The stability equations are 
solved for the spatial amplification rate of each 
unstable frequency. 

(iii) An amplitude ratio for each frequency is then 
calculated by integrating the spatial amplifi- 
cation rate in the streamwise direction on the 
body, i.e. 

to (£) = <i) 

(iv) The n factor is then determined by taking the 
maximum of the above quantity at each stream- 
wise location. 

The major problem with the e n method is that 
the n factor does not represent the amplitude of a 
disturbance in the boundary layer but rather an am- 
plification factor from an unknown amplitude A Q , 
The amplitude A 0 represents the amplitude of a dis- 
turbance of specified frequency at its neutral stabil- 
ity point. Its value is related to the external dis- 
turbance environment through some generally un- 
known receptivity process. As a consequence, the 
value of n which determines transition onset must 
be correlated to available experimental data. Un- 
fortunately, the method suffers from the fact that u 
is not constant; moreover, the method is not useful 
in predicting transition onset on three dimensional 
configurations especially when the crossflow instabil- 
ity is important. 1, 2 Additionally, the e n method re- 
quires the use of several computational tools such as 
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a boundary layer or Navier-Stokes flow solver to cal- 
culate the mean flow and the linear stability solver 
to determine the amplification rate i. The process 
requires significant input from the de signer and gen- 
erally requires substantial knowledge of linear sta- 
bility concepts. Methods based on the parabolized 
stability equations 3 (PSE) are being used to deter- 
mine transition onset but they have not received the 
wide acceptance enjoyed by the e n method. Methods 
based on the PSE also require pre-calculation of the 
mean flow and specification of initial conditions such 
as frequency and disturbance eigenfunctions. Meth- 
ods based on linear stability theory only provide an 
estimation of the location of transition and can pro- 
vide no information about the subsequent turbulent 
flow. 

In this work, a different approach has been devel- 
oped which does not require pre-calculation of the 
mean flow and includes the laminar, transitional, 
and turbulent regions in a single computation. The 
approach employs a two-equation model similar to 
that employed in turbulent calculations. It is based 
on the premise that, if a flow quantity can be writ- 
ten as the sum of a mean and a fluctuating quantity, 
then the exact equations that govern the fluctua- 
tions and their averages are identical irrespective of 
the nature of the oscillations, i.e., laminar, transi- 
tional or turbulent. Moreover, if it is possible to 
model the equations governing the mean energy of 
the fluctuations and their rate of decay (or other 
equations) in such a way that one does not appeal to 
their nature, then the resulting model equations will 
be formally identical. However, the parameters that 
appear in the modeled equations will depend on the 
nature of the fluctuations. As an illustration, let us 
assume that we employ a Boussinesq approximation 
to model the stresses resulting from the fluctuations, 
i.e. 


Tij = -pu'^ = 



2 dU„ 
z ij dx n 



( 2 ) 


where 



dUj | dUj 
dxj 


and turbulence fluctuations but the expressions for 
Pt are quite different because the physics governing 
them is different. 

The goal stated above has recently been 
accomplished. 4 In that work, a new turbulence 
model was developed based on the exact equations 
governing the variance of velocity (kinetic energy) 
and the variance of vorticity (enstrophy). The form 
of the resulting equations was arrived at without 
making use of the fact that the flow was turbulent. 
Therefore, the modeled equations apply to flows of 
any nature, whether laminar, transitional or turbu- 
lent. In order to apply the modeled equations it is 
necessary to specify a “stress-strain relation.” In 
Ref. [4] the set of equations was applied to turbu- 
lent flows. This allowed the model constants to be 
determined by comparing with experimental data. 

The approach followed in this investigation is 
based on the modeled equations of Ref. [4] but with 
different “stress-strain relations.” Further, it is as- 
sumed that the model constants derived in Ref. [4] 
remain unchanged. 

To show the nature of the new “stress-strain rela- 
tions,” Eqs. (25) and (27) of Ref. [4] are re-written 
as 
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p is the density, Ui is the mean velocity, &ij is the 
Kronecker delta and pt is the coefficient of viscosity 
brought about by the presence of fluctuations. The 
form indicated in Eq. (2) is used for both laminar 


Rt — C — ““ ftt — e UmdlU m 

and, 
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(7) 



Lev* of 
Modeling 


UJV 

u! 


3.2 ReJ. 3/2 


where U e is the velocity at the edge of the boundary 
layer, Re$* is the edge Reynolds number based on 
displacement thickness 6 * , and u is the frequency. 

The eddy viscosity resulting from fluctuations in 
the laminar region can be modeled by 


vt — k } Cp — 0.09 (8) 

where r^ t is a viscosity time scale. Using the fre- 
quency of the dominant T-S disturbance, the viscos- 
ity time scale can be modeled in the laminar region 
as 


Figure 1: Schematic of modeling objectives in 

present work 


Ut= min(/* f ,l) (5) 



Tjc is a representative decay time for the kinetic en- 
ergy and ut is the kinematic eddy viscosity. The 
closure coefficients for the model are given in Ta- 
ble 1. 

As may be seen from the governing equations, one 
needs to specify u t and t* to effect closure. Within 
the laminar region, these quantities will be deter- 
mined based on results of linear stability theory. The 
shaded areas of Figure 1 illustrate the objectives of 
the current approach. The diagram shows that the 
exact equations govern the fluctuations regardless of 
their nature, i.e. laminar, transitional, or turbulent. 
The stresses are then modeled by the Boussinesq ap- 
proximation for fill types of fluctuations. The work 
of Robinson and Hassan 5 developed expressions for 
v t and Tfc in the turbulent region. The current work 
develops expressions for v t and r* in the laminar and 
transitional regions. 

For subsonic Mach numbers and regions where 
crossflow instability is unimportant, the dominant 
mode of instability is the first mode or the Tollmien- 
Schlichting (T-S) mode. For low speed flows, the 
dominant disturbance frequency at breakdown is 
well predicted by the frequency of the first mode 
disturbance having the maximum amplification rate. 
Using the work of Obremski et al., 6 Walker 7 showed 
that this frequency can be correlated by 


a 

where a is a model constant, 
region, the representative decay 
energy is modeled as 


(9) 

Within the l amin ar 
time for the kinetic 


r k t 



(10) 


Following the work of Robinson and Hassan 5 , the 
turbulent region time scales are given as 


k 

T » ~ 

(ii) 

VI* 

n 

(12) 


The two regions can be combined by using the 
concept of the intermittency. The intermittency, T 
represents the fraction of time that the flow is tur- 
bulent. At a given point, the flow is laminar (1 - T) 
of the time and turbulent F of the time. This allows 
the viscosity time scale to be written as a transitional 
viscosity time scale, i.e. 


T/i — (1 F)r^ + Fr^ t (13) 

Combining Eqs. 8, 9, and 11a transitional eddy 
viscosity can be written as 

vt = Cufpk (14) 

where 

v = (i - r) ( j) + r (i) (is) 

and 

U = i + r (/„, - 1) (16) 

Following a similar approach, the transitional decay 
time for the kinetic energy can be written 
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- r > (v s ) + r (t) <17) 

The inter mi ttency, T, is currently given by the 
Narasimha et al. 8 expression 

r(i) = l-exp(-.^ 2 ) ( 18 ) 

with 

£ = max(x - x*,0)/A, .4 = 0.412 (19) 

where x* denotes the transition point and A charac- 
terizes the extent of transition. For attached flows, 
an experimental correlation between A and x t is 

Re x = 9.0 He* 75 (20) 

The transition point xt is determined as a part of 
the solution procedure and will be discussed along 
with the results. 

Results and Discussion 

I. Implementation in Boundary Layer Solver 

The present model was initially incorporated into 
the boundary layer code of Harris et al. 9 For such a 
calculation, initial profiles are needed to begin the 
marching procedure. For an initial station s along 
the surface, the dominant disturbance frequency is 
given by Eq. 7. This frequency was then employed 
in the linear stability code of Macaraeg et al. 10 to 
calculate the eigenvalues and eigenfunctions which 
make up the velocity disturbance. The velocity dis- 
turbance corresponding to this dominant frequency 
can written as 

u\ = U{(y) exp [i(as - u it)] + c.c. (21) 

where a, the complex spatial amplification rate, and 
Uj(y), the eigenfunction, are determined from the 
stability code for the specified frequency. The initial 
profile of k is then calculated from 

k = 1«X (22) 

The amplitude of the disturbance was set from the 
specified freestream intensity T u defined as 

<23 > 

Note that in the laminar region, Eq. 3 does not de- 
pend on £. An initial profile for ( is needed for the 


marching procedure but results were found to be in- 
dependent of the choice of the C profile. 

The model constant a is determined in the current 
work by comparing with the flat-plate experiments 
of Schubauer and Klebanoff 11 and Schubauer and 
Skramstad. 12 These classical experiments are well 
documented in the literature and cover a range of 
freestream turbulence intensities. The constant is 
correlated as a function of the freestream intensity 
as 

a = 0.0095 - 0.019 Tu + 0.069(Tu) 2 (24) 

Figures 2 and 3 show the effect of the initial pro- 
files on the onset of transition by plotting skin fric- 
tion as a function of position along the plate. X Q 
is the location where the solution was started and 
X t is the location of minimum skin friction. Fig- 
ure 2 compares the present approach with the ex- 
periment of Ref. [11], while Figure 3 compares with 
the Tu = 0.2 experiment of Ref. [12]. Both figures 
clearly illustrate that the calculated X t values axe 
almost independent of X 0 and are well within the 
scatter of the experimental results. 

The rms amplitude ratio of the velocity distur- 
bance can be obtained from the kinetic energy of 
the fluctuations as 

rr I U a + V * 2 + W r2 

A = vk = y (25) 

For the natural transition process there is a region 
of linear growth, followed by a region of non-linear 
growth which occurs as the amplitude of the velocity 
disturbance becomes sufficiently large. Transition 
occurs after the onset of the non-linear growth and 
this non-linear growth continues through the tran- 
sitional region. After the non-linear growth region 
the modes which make up the disturbance become 
saturated. This saturation of modes characterizes 
the turbulent region. Figure 4 is a plot of the ampli- 
tude ratio as calculated by the present method for 
laminar flow (i.e. T = 0), with k Q representing the 
initial amplitude of k. As can be seen from the fig- 
ure, the present method does predict the expected 
linear, non-linear, and saturated regions. 

The location of minimum skin friction is com- 
monly taken as the onset of transition. In prac- 
tice, it is very difficult to determine the minimum 
skin friction point in an evolving calculation. This 
can be due to either the transient nature of Navier- 
Stokes calculations or due to local oscillations in the 
skin friction itself as seen in later airfoil results. To 
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alleviate this problem, an alternate criteria was de- 
veloped by comparing with the flat plate results. It 
was observed that the present method predicted the 
local skin friction minimum at a location along the 
plate where the maximum -J* « 9%. This criteria 
can be stated differently by noting that the “turbu- 
lent” Reynolds number, Rt, can be written in the 
laminar region (T = 0) as 

Rt = ~ (26) 

Written in this form, Rt can be considered a fluc- 
tuation Reynolds number instead of a turbulent 
Reynolds number. The location of transition on- 
set, xt , is then determined as the minimum distance 
along the surface for which Rt > 1. 

Using this criteria Figure 5 compares the results 
for the present method with the skin friction experi- 
mental data of Ref. [11]. As seen from the figure, the 
present method does a good job of predicting the 
transition onset as well as the transitional and tur- 
bulent regions. Results for the present method are 
summarized in Table 2. As seen from the data, the 
onset of transition predicted by the current method 
compares well with the experimentally observed lo- 
cations. 

To determine the validity of the correlation for the 
model constant, a, the airfoil experiments of Mateer 
et al. 13 were considered. The two-dimensional airfoil 
shown in Figure 6 was used in the experiments over a 
range of angle of attack and Reynolds number. The 
angle of attack considered was -0.5° at Reynolds 
numbers of 0.6, 2, and 6 million. The Mach number 
was 0.2 and the freestream rms pressure and velocity 
fluctuation levels were 0.02goo and O.OOSf/oo respec- 
tively. 

Figures 7 through 9 compare the skin friction re- 
sults for the airfoil with the experimental data and 
Navier-Stokes predictions of Ref. [13]. The Navier- 
Stokes results of Ref. [13] were obtained using the 
e n method and the turbulence model of Menter. 14 
The boundary layer results for the present method 
were computed by using the boundary layer code of 
Harris et al. 9 with a pressure distribution obtained 
from an Euler solver. The initial profile of k for these 
airfoil cases was chosen as the freestream value. 

Figure 7 compares the results of the present 
method for the Re c = 6xl0 5 6 case. Better agreement 
with the experimental data is seen with the transi- 
tion onset prediction of the current method when 
compared with the prediction of the e n method. 
For the airfoil results presented, the e n method did 
not predict transition for the upper surface at all. 
Transition for the e n method was determined on 


the upper surface by the location where the laminar 
boundary layer began to separate. Figure 8 presents 
the results for the Re c — 2xl0 6 case. Again, the 
present method does a better job of predicting tran- 
sition onset. Figure 9 presents the results for the 
Re c — 0.6xl0 6 case. As seen from the data and 
discussed in Ref. [13], the boundary layer is very 
close to separation. Computations with the bound- 
ary layer code are not able to continue past the sep- 
aration point. 

II. Implementation in Navier-Stokes Solver 

The present method has also been incorporated 
in a Navier-Stokes solver. This eliminates the need 
for obtaining the pressure distribution required by 
boundary layer type methods. Additionally, since 
the Navier-Stokes approach is a time marching 
scheme instead of a spatial marching scheme, speci- 
fying initial spatial profiles of k and £ is not possible. 
By using the Navier-Stokes approach, the need for 
linear stability theory to provide an initial profile is 
eliminated. Additionally, it is no longer necessary to 
specify a pressure distribution. 

Figure 5 again compares the results for the present 
method with the skin friction experimental data of 
Ref. [11]. As seen from the figure, the Navier-Stokes 
and boundary layer approaches predict the same 
transition onset location but the Navier-Stokes ap- 
proach predicts a slightly higher value of the skin 
friction in the transition region. 

Figures 7 through 9 also compare the skin friction 
calculated by the present method with the Navier- 
Stokes approach with experiment and other com- 
putations. Figure 7 compares the results for the 
Re c = 6x10° case. The present Navier-Stokes ap- 
proach predicts results nearly identical to the bound- 
ary layer approach. The peak in skin friction on the 
lower surface is slightly higher than the boundary 
layer approach but both methods predict the onset 
of transition much better than the e n method. The 
skin friction in the turbulent region is slightly better 
predicted than the boundary layer approach for the 
upper surface and slightly worse than the bound- 
ary layer approach for the lower surface. Figure 8 
presents the results for the Re c = 2xl0 6 case. Again, 
the present Navier-Stokes approach predicts nearly 
identical transition onset locations when compared 
with the boundary layer approach. The Navier- 
Stokes approach does a slightly better job on the 
upper surface but both Navier-Stokes and bound- 
ary layer approaches over-predict the skin friction for 
the lower surface in the turbulent region. Figure 9 
presents the results for the Re c = 0.6x10 s case. For 
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this case the boundary layer approach was not able 
to proceed past the sepaxation point. The Navier- 
Stokes approach, however, was able to calculate the 
transitional and turbulent regions. In contrast to 
the e n method, the present method is able to predict 
transition onset for the upper surface without fixing 
the location at the separation point. This is because 
the present Navier-Stokes calculations do not predict 
separation on the upper surface. Transition onset in 
the present method is fixed by the separation point 
for the lower surface. 

Conclusions 

The current results demonstrate that the present 
approach addresses those aspects of the breakdown 
process that follow the growth of linear disturbances. 
The e n method is based on the conjecture that there 
is a critical amplification factor of the Tollmien- 
Schlichting waves at transition. The current ap- 
proach does not support this conjecture. A key 
finding of the current results is that the most am- 
plified frequency plays a major role in determining 
the onset of transition. Current results show that 
the method is an inexpensive alternative to the e n 
method and provides promising transition onset pre- 
diction capabilities for a range of Reynolds numbers 
and freestream intensities. 
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Tables 


Constants 

fc - < ii 

C u 

0.09 1 

K 


at 3 

0.35 l| 

A 

HUH 

A 

■zai 

A 

2.37 

IKM 


A 

0.75 | 

IHK9HH 

msm 

CTr 

2.00 

X 

1.80 

** 

1.46 


4.0 


msm 

* i 

0.1 || 


Tkble 1: A; — £ Model Constants 


Case 

Tu 


X t (Current Method) 

Error (%) 

Schubauer-Klebanoff 

0.03 

5.26 ft. 

5.39 ft. 


Schubauer-Skramstad 

0.042 

5.25 ft. 

5.24 ft. 


Schu b auer- Skrams tad 

■ 

5.08 ft. 

4.86 ft. 


Schubauer-Skramstad 

■ 


4.08 ft. 

0.74 

Schubauer-Skramstad 


3.32 ft. 

3.41 ft. 


Schu bauer-Skramstad 

0.34 

2.58 ft. 

2.52 ft. 



Tkble 2: Flat Plate Experiments, M re f = 0.071, 
Pref = 1.01325x10 s Pa, T ref = 293 K. 
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X(ft.) 

Figure 2: Effect of Initial Solution on Transition 
Onset, Comparison with Experiment of Ref. [11] 



Figure 4: Amplitude ratio versus distance along 
plate. Schubauer-Klebanoff Experiment 



Figure 3: Effect of Initial Solution on Transition 
Onset, Comparison with Experiment of Ref. [12] 



Figure 5: Comparison of present method with 

the experiment of Schubauer and Klebanoff, Re = 
1.67x10* /m, Navier-Stokes &c Boundary Layer 
Codes 
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Figure 6: Airfoil geometry for experiment of Ma- Figure 8: Comparison of present method and e n 
teer et al. i3 method with the airfoil experiment of Mateer et al. i3 

, Re c = 2 million 
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Figure 7: Comparison of present method and e n 
method with the airfoil experiment of Mateer et al. 13 
, Re e = 6 million 



04 QlS 04 Out 04 14 


Figure 9: Comparison of present method and e n 
method with the airfoil experiment of Mateer et al. 13 
, Re c =s .d million 
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Abstract 

An approach developed earlier by the authors for 
determining transition onset resulting from instability 
of Tollmien-Schlichting waves is extended to crossflow 
instabilities. In this approach, a two-equation model 
similar to that used in turbulent flows is employed. 
The theory takes into consideration the role of two 
major environmental influences: turbulence intensity 
and surface roughness. Comparisons with data for an 
infinite swept wing and an infinite swept plate under 
the influence of a favorable pressure gradient shows 
excellent agreement with experiment. Calculations 
employed a boundary layer code: as a result, compu- 
tational costs axe minimal. 

Introduction 

I T is well known 1,2 that the dominant mode of in- 
stability for the majority of swept-wing flows is the 
crossflow instability. The boundary layer on wings 
of moderate or greater sweep generally contains sig- 
nificant crossflow. The velocity profile in this case 
can be separated into a component in the stream- 
wise direction and a component in the crossflow di- 
rection. Since the crossflow profile always contains an 
inflection point, a strong inflectional instability is ex- 
pected in regions where the crossflow velocity increases 
rapidly. 3 This increase occurs in regions of negative 
pressure gradient, e.g. near the leading edge. In 
this favorable pressure gradient region, streamwise or 
Tollmien-Schlichting (T-S) instabilities are stabilized 
and the crossflow (CF) instability dominates the tran- 
sition process. 

Methods employed in transition prediction are dis- 
cussed in detail by Haynes et al. 4 The practical meth- 
ods that are in current use are the e n method, which 
is based on linear stability theory (LST) and methods 
based on the parabolized stability equations (PSE). 
It is generally concluded that the e n method is not 
suited for predicting transition onset resulting from 
CF instabilities. This is because these instabilities 
are sensitive to surface roughness, and because they 
are dominated during their development by large non- 
linear effects that extend over a wide region. As a 
result, growth rates are not well predicted. Even if 
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such rates were known, the variety of methods used to 
calculate the n factor result in a great deal of scatter. 3 
The nonlinear PSE method 5 does a much better job of 
predicting growth rates if the disturbance inputs are 
known. Because such inputs for flight conditions and 
wind tunnels are not well known, its success in pre- 
dicting transition onset from CF instabilities is yet to 
be demonstrated. 

In most cases, the only data available is the turbu- 
lence intensity and very little information is provided 
on amplitude or frequencies of the various modes that 
may be present. In addition, stability analyses make 
use of a number of codes and require a great deal of 
knowledge and skill on the part of the user. Rec- 
ognizing the importance of determining the onset of 
transition over a wide range of operating conditions 
for a vehicle designer, and in an effort to simplify 
the procedure for calculating such onset, an approach 
similar to that used in studying turbulence was de- 
veloped by Warren and Hassan. 6 It was demonstrated 
in situations where T-S waves played the major role 
in affecting transition. The approach required knowl- 
edge of turbulence intensity and was implemented in 
both boundary layer and Navier-Stokes codes. Its 
development was based on the observation that, if 
a flow property can be written as the sum of mean 
and fluctuating components, then the exact equations 
governing the variance of velocity (kinetic energy of 
fluctuations) and the variance of vorticity (enstrophy) 
are the same irrespective of the nature of the fluctu- 
ations. Moreover, if such equations can be modeled 
without declaring the nature of fluctuations, the re- 
sulting equations are capable of describing all types of 
fluctuations in a fluid flow, i.e. laminar, transitioned 
or turbulent. What distinguishes one type of fluctua- 
tion from another is the appropriate stress-strain law 
required to close the resulting system of equations. 

The objective of this work is to extend the approach 
of Ref. 6 to crossflow instabilities. In order to incorpo- 
rate the procedure within existing computational fluid 
dynamics (CFD) codes, we opted to employ an eddy 
viscosity model. For such models, eddy viscosity can 
be written as 

v t oc k Tt r (1) 

where k is the fluctuation kinetic energy per unit mass 
and r tr is an appropriate time scale. In Ref. 6, r tr 
was related to the frequency of the most amplified 
waves. In the present case, r tr was obtained from the 


1 OF 9 


American Institute of Aeronautics and Astronautics Paper 97-2245 



observation that the wavelength of stationary cross- 
flow vortices lies in the range of three to four times 
the boundary layer thickness, 6. As a result, r tr was 
chosen as 

3.5J 

Ttr * ~qT ( 2 ) 

where Q e is the streamwise velocity at the edge of the 
boundary layer. As is shown below, the proportion- 
ality model constant depends on both the turbulent 
intensity and surface roughness. These are the two 
most dominant environmental factors that affect tran- 
sition. 1 

Determining the onset of transition does not pose 
a problem from a computational standpoint; it does, 
however, represent a problem from an experimental 
standpoint. The data used in the validation of this 
method is that of Saric 1 and his colleagues and of 
Bippes 1 and his colleagues. The first group employed 
napthalene flow visualization while the second group 
defined onset where the intermittency is 50%. Because 
visualization data is somewhat subjective, the latter 
criterion is used in this work to facilitate comparison 
with experiment. 

There is one important difference in implementa- 
tion between this work and that of Ref. 6. In Ref. 6, 
the method was implemented in both boundary layer 
and Navier-Stokes codes because skin friction data 8 
was available to validate both transition onset and 
flowfield calculations. Thus, the codes calculated the 
laminar, transitional, and turbulent regions in a sin- 
gle computation without interference from the user. 
Because flows being considered were attached, imple- 
mentation was limited to a three-dimensional bound- 
ary layer code suited for calculation of infinite wings. 
The code is a modification of an earlier code by Harris 
and Blanchard. 9 

Model Formulation 

The present method takes advantage of the fact that 
whether one deals with transition or turbulence, each 
flow quantity is set as a sum of a mean and fluctuating 
component. Moreover, the exact equations that gov- 
ern the fluctuations are the same irrespective of the 
nature of the fluctuations. In the work of Robinson 
and Hassan, 10 the exact equations governing the vari- 
ance of velocity (kinetic energy of fluctuations) and the 
variance of vorticity (enstrophy) were modeled with- 
out making use of the fact that the flow was turbulent. 
Therefore, the modeled equations apply to flows of any 
nature, whether laminar, transitional, or turbulent. In 
order to apply the modeled equations, it is necessary 
to specify a "stress-strain relationship.” Following the 
same method used in Warren and Hassan, 6 the present 
approach for including crossflow instabilities is based 
on the modeled equations of Ref. 11 with different 
stress-strain relations. ? The equations for the kinetic 


energy and enstrophy of Ref. 11 are re-written for low- 
speed and favorable pressure gradient flows as 
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t* is a representative decay time for the kinetic energy 
and v t is the kinematic eddy viscosity. The closure 
coefficients for the model are given in Table 1. 

As may be seen from the governing equations, one 
needs to specify i/ t and r* to effect closure. Following 
the same approach as Ref. 6, the current work develops 
expressions for v t and r* in the laminar region. Using 
the concept of intermittency, the laminar expressions 
for v t and r * are blended with the turbulent expres- 
sions developed in the work of Robinson and Hassan. 12 
A schematic of the approach is shown in Figure 1. 

The eddy viscosity resulting from fluctuations in the 
laminar region is written as 


Vt^C^k r U( , = 0.09 (6) 

where r Ut is a viscosity time scale which characterizes 
the dominant disturbance in the laminar region. Ref. 6 
only considered the streamwise or T-S instability. Us- 
ing a correlation by Walker 13 for the dominant T-S 
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( 10 ) 


Constants j A: — C | 

I 0.09 I 

K 

0.40 | 

<*3 

0.35 

04 

0.42 

$5 

2.37 

06 

0.10 

$7 

1.50 

8$ 

2.30 


0.07 

I 

1.80 

J_ 

*C 1 

1.46 

6 

0.1 


Table 1 k — £ model constants 


/ , r\ 'M Jr 

r MrcF — ( a + /) g 

where a is the model constant given by Eq. (8) and / is 
a model constant which depends upon surface rough- 
ness. Q e is the resultant boundary layer edge velocity. 
A cf is the wavelength of the crossflow disturbances 
and has been shown by Arnal 2 and other researchers 
to be between three and four times the boundary layer 
thickness, 5. Additionally, the work of Muller and 
Bippes 15 demonstrated that the wavelengths of both 
traveling and stationary crossflow disturbances were 
the same but shifted in phase by Acf/2. Choosing 

A cf ~ 3.5 <5 

the viscosity time scale in the laminar region is mod- 
eled in the present work as 



A 

j Lftvat of 
Modeling 


Fig. 1 Schematic of modeling objectives in present 
work 

disturbance frequency, tu, the dominant time scale for 
the T — S disturbances was modeled as 

r L*e,TS ~ ~~ (7) 

tJ 

where a is a model constant given in Ref. 6. The ex- 
pression for a is rewritten here as 


„3.55 

r M.CF "" “b /) Q 

The expression for 7> t is the same form as in Ref. 6 
with the addition of the roughness constant /, 


— = (a + /)-S (12) 

T k e V 

Following the work of Robinson and Hassan, 11 the tur- 
bulent region time scales are given as 



(13) 


As in Ref. 6, the laminar and turbulent regions are 
combined through the intermittency, T. This allows a 
transitional eddy viscosity to be written as, 


i/t — Qn k Tjj 


(14) 


where 


= (1 - n 


(a + /) 


3.5 S' 
Qe . 



(15) 


a = 0.00819 + 0.069 (Tu - 0. 138) 2 (8) 

and the freestream intensity, Tu, is defined as 


Tu = 100 



(9) 


In contrast to Ref. C, the current work develops 
a viscosity time scale which characterizes only cross- 
How disturbances. As discussed in Refs. 1, 2, 7, and 
14, the wavelength of the dominant crossflow distur- 
bance varies with the boundary layer thickness, 6. In 
addition to freestream disturbance levels, crossflow in- 
stabilities are sensitive to surface roughness. Using 
this information, the present method models the vis- 
cosity time scale for the crossflow instabilities as 


Following a similar approach, the transitional decay 
time for the kinetic energy can be written 

i = (i-n[<. + /)£s] + r(£) (16) 

The intermittency, T, is currently given by the 
Narasimha et al. 16 expression 

r(z) = 1 -exp(-0.412£ 2 ) (17) 

with 

£ = max(x — x t > 0)/A (18) 

where A characterizes the extent of the transitional 
region. An experimental correlation between A and x t 

is 
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Re x =9.0Re°™ (19) 

x t in the above expression represents the location 
where turbulent spots first appear. The present work 
determines x t by the same method used in Ref. 6. In 
that work, it was shown that the minimum skin fric- 
tion occurred when the fluctuation Reynolds number, 
Rt > exceeded one. Rt is defined as 


Thus x t is the minimum location where Rt > 1. 

Results and Discussion 

The present method is compared with the swept- 
wing experiments of Dagenhart et al. 14 and Radeztsky 
et al. as well as the swept flat plate experiments of 
Muller and Bippes 15 and Deyhle and Bippes. 7 Each 
of these experiments simulated infinite wing flow con- 
ditions. The infinite swept wing assumption of zero 
spanwise gradients was applied to the governing three- 
dimensional boundary-layer equations 18 and the equa- 
tions governing the fluctuations, Eqs. (3) and (5). The 
resulting equations were then solved using a modified 
version of the boundary-layer code of Harris and Blan- 
chard. The necessary pressure distribution was given 
in the swept plate experiments and was computed from 
an Euler solver for the swept wing cases. The initial k 
profile was set from the freestream intensity, Tu, of the 
experiments. 6 Figure 2 shows the coordinate system 
used for the swept wing and swept plate geometries. 

The infinite swept-wing experiments of Dagenhart et 
al. and Radeztsky et al. 17 were carried out on a 45° 
swept wing with a NLF(2)-0415 airfoil cross-section at 
a 4° angle of attack. Figure 3 illustrates the NLF(2)~ 
0415 profile and compares the experimental pressure 
coefficient on the upper surface with the computed re- 
sults used in the present study. Transition onset was 
determined by napthalene flow visualization. As de- 
scribed in Dagenhart. 19 transition was determined by 


observing the sublimation pattern of the napthalene 
spray on the painted airfoil surface. The sublima- 
tion indicated higher shear levels characteristic of the 
transitional and turbulent regions. Traditionally on- 
set of transition is defined by a minimum skin Auction 
coefficient or minimum heat transfer coefficient. How- 
ever, for comparisons with available data, it is assumed 
that the experimentally observed transition points re- 
ported in Ref. 17 actually occurred at locations wit hin 
the transitional region, i.e. non-zero intermittency. 
The swept flat plate experiments of Refs. 7 and 15 
measured intermittency and defined transition as the 
location where T = 0.5. It is assumed, for convenience, 
that the data reported in the work of Radeztsky et al. 17 
corresponded to a location where T = 0.5 as well. 

The experiments reported in Radeztsky et al. 17 were 
carried out over a range of Reynolds numbers with 
three different model surface finishes. The 9/rm case 
was the original painted finish of Dagenhart et al. 14 
and is a “peak-to-peak” measured roughness of the 
surface finish. 17 The 0.25pm and 0.5pm cases are root- 
mean-squared (rms) measured roughness of the surface 
finish. By comparing with these cases for a given chord 
Reynolds number, the model constant / in Eq. (11) is 
correlated with roughness as is given as 


/ = 0.003 



Ren - 0.92 


where 


( 21 ) 


h ref = 1 pm 


t'oo 

h is the "‘peak-to-peak” distributed roughness level 
(finish). For rms measured levels, a sinusoidal dis- 
tribution is assumed and the level is increased by a 
factor of \/2. These calculations assume a freestream 
intensity level, Tu = 0.09, reported by Dagenhart. 19 

Figure 4 compares the experimental transition loca- 
tions with the location of T = 0.5 as computed by the 
present method. Excellent agreement is observed for 
all Reynolds numbers and for all three surface finishes. 

To further validate the expression for the model con- 
stant /, the present approach was used to compare 
with the swept flat-plate experiments of Muller and 
Bippes, 15 and Deyhle and Bippes. 7 The reported ex- 
periments were carried out in three different wind tun- 
nels. designated as NWB, 1MK, and NWG. The NWB 
had a freestream intensity level of Tu = 0.08, the 1MK 
had values of Tu — 0.15 and Tu = 0.27 depending on 
whether or not a screen was present, and the NWG had 
an intensity of Tu — 0.57. Several plates of varying 
surface roughness were tested. Transition was deter- 
mined in the experiment by increasing the effective 
freestream velocity until an intermittency of 50% was 
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Facility 

Tu 

Plate Surface 

Experimental 
Rer.tr (xlO 5 * ) 

Present 
Re T ,tr (xlO 5 ) 

Error (%) 

NWB 

0.08 

Wooden plate, R z = 6^m 

6.5 

8.04 

23.7 

1MK 

0.15 

Plate with sandpaper. R z = 40^m 

6.8 

6.12 




Wooden plate, R z = 6/im 

7.5 

7.49 

0.13 



Aluminum plate, sanded. R z = b^irn 

7.7 

7.63 

0.9 



Aluminum plate, polished, R z = 1.8 /im 

8.3 

8.48 

2.17 

lMK/screen 

0.27 

Aluminum plate, polished, R z = l.&fim 

7.8 

7.11 

8.85 

NWG 

0.57 

Aluminum plate, sanded, R z = 5 \im 

5.4 

2.18 

59.6 


Table 2 Comparison of Re Xcttr values at 50% intermittency predicted by the present method and measured 
in the swept plate experiments of Deyhle and Bippes. 7 


observed at the measurement station x c /c= 0.9. The 
crossflow instability was isolated in the experiment by 
imposing a negative pressure gradient along the plate. 
The pressure coefficient measured in the experiment is 
shown in Figure 5 for the three tunnels. 

Table 2 compares the predictions of the current 
method with the experimented results given in Ref. 7. 
As seen from the table, the present approach does an 
excellent job of predicting the transition locations for 
the 1MK tunnel. Transition locations are predicted 
within experimented error for ail intensities and sur- 
face conditions in the 1MK tunnel. The predictions 
by the present method were not as good for the NWB 
case and the transition location for the high-intensity 
NWG case was severely underpredicted. As discussed 
in Ref. 7, the Tu = 0.57 case could have been char- 
acterized by by-pass effects. It is also possible that 
the correlation for a, which was was obtained based 
on transition resulting from T-S instabilities, may not 
be appropriate for CF instabilities. We have no logical 
explanation as to the cause of discrepancy in the NWB 
tunnel. It is worth mentioning that the environmental 
conditions in the tunnel are comparable to those of the 
Arizona State University tunnel. 14,17 

Figure 6 compares the calculation of the intermit- 
tency, T, predicted by the present method at x c /c=0.9 
versus Reynolds number. As seen from the figure, the 
present method compares well with the experimental 
data. 

In Ref. 15, first turbulent spots were detected at 
location x c /c=0.95 for the 6/im plate in the 1MK 
tunnel at Tu — 0.15 and Q ^ = 19m/s. The loca- 
tion of first turbulent spots is the definition of x t in 
Eq. (18) and the present method predicts this location 
at x c /c=0.949 giving excellent agreement with Ref. 15. 

The rms amplitude of the total velocity disturbance 
can be obtained from the kinetic energy of the fluctu- 
ations as, 

t m u t2 + v i2 +w r2 

,4 = Vk = y (22) 

Figure 7 is a plot of the amplitude ratio for the 
Tu = 0.15 and 6 /im roughness case in the 1MK 


tunnel. k 0 in the figure represents the maximum ki- 
netic energy of the fluctuations at the initial station 
in the boundary-layer marching procedure (freestream 
k). As seen from the figure, the present method pre- 
dicts a region of linear growth and a significant region 
on non-linear growth as the onset of transition is ap- 
proached. The disturbance level tends to saturate as 
onset is approached but resumes exponential growth 
when the transitional region is entered. This agrees 
qualitatively with the large eddy simulation (LES) re- 
sults presented by Huai et al. 20 The ability to capture 
the nonlinear growth shown in Figure 7 demonstrates 
the sound basis of the present method. 

Conclusions 

It is shown in this investigation that treating tran- 
sitional flows in a manner similar to turbulent flows 
represents a convenient and an inexpensive method for 
determining transition onset. The key to this approach 
is to determine the relevant eddy viscosity suited for 
describing the underlying physics of the problem. 

Because of the turbulence-like approach developed 
here, the manner in which the various modes interact 
with each other and with the environment to bring 
about transition is integrated out of the problem. Ev- 
idently, one of the significant outcomes of the inter- 
action is the fact that the CF disturbances have a 
wave length that varies between three to four times 
the boundary layer thickness. 
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Fig. 3 NLF(2)-0415 airfoil geometry and comparison of experimental and computed pressure coefficient 

on the upper surface. 
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Fig. 4 Comparison of the present method with the experimental data of Radeztsky et al. 17 
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Fig. 5 Experimental pressure coefficient along the swept plate of Deyhle and Bippes 7 
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Fig. 6 Comparison of the intermittency predicted by the present method with the experimental data of 
Deyhle and Bippes, 7 z c /c=0.9. 
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ABSTRACT 

A unified approach which makes it possible to deter- 
mine the extent and onset of transition in one calculation 
is presented. It treats the laminar fluctuations in a man- 
ner similar to that used in describing turbulence. As a 
result, the complete flowfield can be calculated using ex- 
isting CFD codes and without the use of stability codes. 
The method is validated by comparing the results for flat 
plates, airfoils, and infinite swept wings with available ex- 
periments. In general, good agreement is indicated. 

INTRODUCTION 

Traditionally, the transition problem has been treated 
as a combination of two problems. The first deals with 
the extent of the transition region given the onset, while 
the second deals with transition onset. One of the methods 
employed in calculating the extent is to set the effective vis- 
cosity, /x, in a boundary layer (BL) code or Navier-Stokes 
(NS) solver as 

M = /x* + I (1) 
where subscripts l and t designate laminar and turbulent 
Sows respectively and T is the intermittency or the fraction 
of time the flow is turbulent at a given location. The most 
widely used expression is that of Dhawan and Narasimha. 1 

There are many ways that are being used to specify 
transition onset: arbitrary selection, experimental correla- 
tions, or use of stability theory. Methods based on stability 
theory employ the e n method or a method based on the 
parabolized stability equations (PSE). Methods based on 
stability theory have shown a great deal of promise when 
stream wise or To llmien -S chlich tmg (T-S) waves are the 
do m i nan t cause of transition. The same cannot be said 
when transition is a result of crossflow (CF) instabilities 
because such instabilities are dominated by nonlinear ef- 
fects and surface roughness. An excellent recent review of 
these methods and their limitations is given by Haynes et 
al. 2 

In situations where transition onset is specified from re- 
sults of an experiment, Eq. (1) does not perform well. One 


of the reasons for this behavior is because the above fc 
mula does not account for the laminar fluctuations th; 
eventually lead to transition. This led Young et al. 3 ar 
Warren et al. 4 to employ an expression for /x given by 

a = t*i + [(i - r )m f + rv*] (: 

where fit/ is the contribution of the laminar fluctuation 
The expression for fitf was correlated by using resul 
from linear stability theory for both low and high spet 
flows. Much better agreement with experiment was inc 
cated when Eq. (2) was employed. 

The inability of stability theory to cope with crossflo 
instabilities led Warren and Hassan 5,6 to develop a ne 
approach for determining transition onset. This approax 
is centered around the determination of /n*/. When \nf 
known, then the onset of transition, which may correspor. 
to minimum skin friction, minimum heat flux, or son 
other criterion specified by the user, can be determined ; 
part of the solution. Such an approach then removes tt 
need for stability codes to predetermine transition onse 
Moreover, it addresses transition onset and extent as or 
and not two separate problems. 

In both transition and turbulence each flow quantity 
set as a sum of a mean and fluctuating quantity. Th 
exact equations governing fluctuations are the same. I 
transition, attention is focused on individual modes an 
frequencies with the growth rates of such modes playin 
a crucial role in determining transition onset. In turbr 
lence, equations governing fluctuations are used to deriv 
equations for the mean energy of the fluctuations and it 
dissipation rate. As a result, individual modes do not pla 
any direct role in turbulence calculations. The resultin 
equations governing turbulence are not closed, thus nece 
sitating assumptions on a stress-strain law. 

The approach presented here takes advantage of similar 
ities between laminar and turbulent fluctuations, i.e. th 
exact equations governing energy and it dissipation rat 
are identical. Moreover, it is possible to model these ex 
act equations 7 without invoking the nature of fluctuations 
In order to close the model equations, it is necessary t 
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provide appropriate stress-strain laws. To facilitate imple- 
mentation in existing CFD codes, an eddv viscosity model 
will be employed. 

In the present work, the stress-strain laws for the lam- 
inar fluctuations are derived from observed or computed 
characteristics of T-S and CF disturbances. Because mech- 
anisms responsible for transition are different for the two 
types of instabilities, corresponding stress-strain laws are 
different. In both cases, however, puf is set as 

Ptf — Cjj. p k t (3) 

where k is the fluctuation kinetic energy per unit mass, p is 
the density, and r is an appropriate time scale characteris- 
tic of the type of instabilities begin considered. Although 
the present approach makes no direct use of stability codes, 
expressions for r were arrived at from consideration of re- 
sults obtained from stability theory or from experimental 
observations. 

It is to be noted that this is not the first attempt at 
using equations similar to the turbulence equations to pre- 
dict transition onset. Typical of these attempts is the work 
of Wilcox. 8 " 10 An earlier attempt 9 tried to take advan- 
tage of the results of linear stability theory. In the latest 
effort, 10 the approach used in Ref. 9 was abandoned in fa- 
vor of an approach in which the model constants in the 
k — w model were replaced by functions of the turbulence 
Reynolds number. The functions were selected in such a 
way so as to reduce to the original model constants at 
high turbulence Reynolds number and to reproduce tran- 
sition onset for an incompressible flow over a flat plate. 
As formulated, the model is insensitive to the modes of 
transition 11 and thus will not yield good results outside 
the range for which it was formulated. 12 This may be con- 
trasted with the present model in which separate equations 
are used to model the “laminar” fluctuations and where 
the stress-strain laws that govern those fluctuations are 
not only sensitive to the various modes of tr ans ition but 
are dictated by them. 

PROBLEM FORMULATION 


The basis for the present approach is the exact equations 
governing k and £, the enstrophy or the variance of vortic- 
ity. The turbulence model employed in this work is 
that of Robinson and Hassan. 13 The equations governing 
k and £ are those given in the work of Robinson et al. 13 
and can be rewritten for low speed flows as 
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and, 


T ij — P u 'i u 'j — Mt ^2 Sij — - dij ^ — —pkSij (6) 


^ is a representative decay time for the kinetic energy and 
i/t is the kinematic eddy viscosity. The closure coefficients 
for the model are given in Table 1. 


Table 1 k — £ model constants 


Constants 

*-c 

c u 

0.09 

K 

0.40 

Ot3 

0.35 

04 

0.42 

0* 

2.37 

\ 0 8 

0.10 

01 

1.50 

08 

2.30 

<?r 

0.07 

- °V 

0.065 

<Tb 

1.80 

JL 

1.46 

1 <5 

0.1 


As may be seen from the governing equations, one needs 
to specify v t and 7 * for each mechanism to effect closure. 
When T-S waves are considered, v t is chosen as 3 


£1 

Dt 


i§W 

sflcf 

tmij 

9 ( U 'm“0 dk 

dxi dx m 

+ a*,[ 

K)£]- 


+ bij + Sij ~ 


08 


** + <T 

pkSl 


(5) 
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where 


r w =r, 
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ui is the frequency of the first mode disturbance having 
the maximum amplification rate and a is a model constant 
that depends on the freestream intensity, Tu, defined as 


Tu = 100 



(9) 
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In the above expression, Qqo is the freestream velocity. 
The frequency u is given by a correlation developed by 
Walker^ 4 as 

^ = 3.2 ReJ. 3 ' 3 (10) 

where Q e is the velocity at the edge of the boundary layer 
and .Re*, is the edge Reynolds number based on a displace- 
ment thickness J*. The model constant a was correlated 
using flat plate experiments by Schubauer and Klebanoff 15 
and Schubauer and Skramstad 16 as 0 

a = 0.069 (Tu - 0.138) 2 + 0.00819 (11) 

Within the laminar region, the representative decay time 
for the kinetic energy is 

— = a—S (12) 

Tk, V 

Similarly, when transition is a result of CF instabilities 6 

^=Wr=(a + /)^ (13) 

where A cf is the wavelength of crossflow disturbances. 
Based on numerous experimental and computational re- 
sults which found wavelengths of CF disturbances between 
3 and 4 times the boundary layer thickness, we use here 

A cf = 3.5 8 (14) 

Since stationary CF disturbances are generated by sur- 

face roughness and traveling disturbances are generated by 
freestream disturbances and surface conditions, a correla- 
tion reflecting the influence of surface conditions must be 
included in the model. Using one of the sets reported by 


Radeztsky et al., 17 / was correlated as 6 


/=o - oo3 [(£) 

-0.8 

Reh - 0.92 

(15) 

where 



Kef = 

lfj.ni 

(16) 

Re h = 

Qooh 

I'oo 

(17) 


and h is the “peak-to-peak” distributed roughness level. 
The decay time for the kinetic energy is chosen as 


~ = (a + f)^S ( 18 ) 

Tk t V 

Following the work of Robinson and Hassan, 13 the tur- 
bulent region time scales are given by 

7^ = r M< = A (19) 

Because the thrust of this work is the prediction of tran- 
sition onset, a simple intermittency approach is used in the 
transition region. As a result 

T u = (l-r)r w +rr M , (20) 

i = (i-r)-J- + r— (2i) 

** Tkt Tkt 


The intermittency, T, is given by Dhawan anc 
Narasimha’s 1 expression, i.e. 

r(x) = 1 - exp (— 0.412£ 2 ) (22) 

with 

f = max(x - x t ,0)/A (23) 

A is a characteristic extent of the transitional region. Ar. 
experimental correlation between A and x t is 

Re x = 9.0 Re° x J 5 (24 

with x t being the location where turbulent spots first ap 
pear, or where skin friction is minimum It is shown in tin 
work of Warren et al. 5 that this location is well representec 
by the relation 

R t = - Li! = 1.0 (25 

Up V 

Thus, x t is the minimum location where Rt > 1. 

RESULTS AND DISCUSSION 

Results presented here employed both boundary layei 
and Navier-Stokes codes. The boundary layer (BL) code 
is an adaptation of the code developed by Harris anc 
Blanchard. 18 In addition to incorporating the present tran 
sition /turbulence model, it was extended to handle infinite 
swept wings. The Navier-Stokes (NS) code is an adapta 
tion of an earlier code developed by Gaffney et al., 19 whicl 
employs an upwind Roe scheme and four step Runge-Kutt: 
time stepping method. When employing a boundary laye 
code, the pressure distribution is obtained from an Eule: 
code or experiment. 



X ( fL ) 

Figure 1 Comparison of present method with the 
experiment of Schubauer and Klebanoff, Rt = 1.67 x 
10 6 /m, Navier-Stokes & Boundary Layer Codes 

Transition resulting from T-S instabilities is discussed 
first. The only environmental parameter that needs to be 
specified is the freestream intensity. Figure 1 compares 
predictions of present theory with measurements over a 
flat plate. 15 With the current model, the only environmen- 
tal condition that needs to be specified is the freestream 
turbulence level. Both computational schemes give similar 
results with regard to onset and extent of transition. The 
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Table 2 Flat Plate Experiments, AL e / = 0.071, P re f = 1 .01325 x 10 5 Pa, T r ef 293 K. 


Case 

Tu 

X t (Experiment) 

X t (Current Method) 

Error (%) 

Schubauer-Klebanoif 

0.03 

5.26 ft. 

5.39 ft. 

2.47 

S chu b auer- Skr ams t ad 

0.042 

5.25 ft. 

5.24 ft. 

0.19 

Schubauer-Skramstad 

0.10 

5.08 ft. 

4.86 ft. 

4.3 

Schubauer-Skramstad 

0.20 

4.05 ft. 

4.08 ft. 

0.74 

Schubauer-Skramstad 

0.26 

3.32 ft. 

3.41 ft. 

2.7 

Schubauer-Skramstad 

0.34 

2.58 ft. 

2.52 ft. 

2.3 


transition onset locations predicted by the present method 
are compared with the flat pate experiments of Schubauer 
and Klebanoff 15 and Schubauer and Skramstad 16 in Ta- 
ble 2. As seen from the table, excellent results are obtained 
for all ff eestream intensities reported. The next set of com- 
parisons involve the data of Mateer et ail. 20 They presented 
skin friction measurements over the supercritical airfoil 
shown in Figure 2 for a freestream Mach number, A/*,, of 



Figure 2 Airfoil geometry for experiment of Mateer 
et al. 30 

0.2 and a range of Reynolds numbers and angles of attack, 
a. Moreover, they compared their measurements with the 
e n method. Emphasis will be placed on comparisons with 
a = -0.5° cases because this is the angle of attack where 
large discrepancies between the e n method and experimen- 
tal data were noted. It is suggested that the e n method 
is incapable of predicting transition for these cases when 
transition is dominated by Reynolds number effects and 
not determined solely by laminar separation. 20 

Figs. 3-5 show a comparison of predictions of present 
theory with experiment and calculations reported in the 
work of Mateer et al. 20 Both BL and NS calculations are 
shown. It is to be noted that when transition is a result 
of flow separation, BL calculations are terminated at sep- 
aration. Fig. 3 compares results for Re c = 0.6 x 10*. For 
this case the e n method sets transition on both upper and 
lower surfaces as the location of laminar flow separation. 
As may be seen from the figure, the present theory gives 
good agreement with experiment for both upper and lower 
surfaces. For Re c = 2 x 10*, Fig. 4 shows that the present 
method gives much better agreement than the e n method 



Figure 3 Comparison of present method and e n 
method with the airfoil experiment of Mateer et al. 30 
, Rec — -6 million, a = —0.5° 



X/C 


Figure 4 Comparison of present method and e n 
method with the airfoil experiment of Mateer et al. 20 
, Rt c = 2 million, a = —0.5° 

with experiment for both upper and low surfaces. A sim- 
ilar conclusion is reached for Re c = 6 x 10* as shown in 
Fig. 5. It is to be noted however, that for this case, the cur- 
rent model overpredicts skin friction over part of the airfoil 
while giving reasonable estimates of onset location. It is 
difficult to pinpoint the cause of the discrepancy in skin 
friction for this case especially when good agreement for 
the other cases was indicated. Possible contributing fac- 
tors may be the expression used to describe intermittency 
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Figure 5 Comparison of present method and e n 
method with the airfoil experiment of Mateer et al. 20 
, Re c = 6 million, a = -0.5° 


NS Results of Msteer. et eJ. 



Figure 6 Comparison of present method and e n 
method with the airfoil experiment of Mateer et al. 20 
, Re c = 2 million, a = 3.5° 


or increased blockage in the t unn el 
For higher angles of attack, the present method and 
e n methods are somewhat comparable in their predic- 
tions. Figure 6 compares calculations for a = 3.5° and 
Re c = 2 x 10 6 . It is to be noted that both models over- 
! estimate skin friction on the upper surface. For the upper 

' surface, both models predict transition at the location of 

\ separation. However, for the lower surface the present 

method predicts a location slightly upstream of the loca- 
tion of laminar separation. For this case, the e n method 
predicts transition at the location of laminar separation for 
both upper and lower surfaces. 

The next set of comparisons address transition resulting 
from crossflow instabilities. Figure 7 shows the coor dina te 
system used for the swept wing and swept plate geome- 
tries. For these cases, available data give tr ans ition onset 
locations . 17 * 21 “ 23 Unfortunately, s kin friction data is not 
provided. Moreover, flow separation does not play any 
role in determining transition onset in available data. As 
a result, computations presented for CF instabilities were 



Figure 7 Coordinate system on the swept wing ai 
flat plate 






Figure 8 NLF(2)— 0415 airfoil geometry and compai 

son of experimental and computed pressure coefficie? 
on the upper surface. 

based on a three-dimensional boundary layer code suite 
for calculating infinite swept wings. For such flows, tl 
boundary layer equations do not depend on the spanwr 
coordinate. 24 

The infinite swept wing experiments of Dagenhart < 
al. 21 and Radeztsky et al. 17 use a 45° swept wing with 
NLF(2)-0415 cross-section at —4° angle of attack. Figure 
illustrates the NLF(2)-0415 profile and compares the & 
perimental pressure coefficient on the upper surface wit 
the computed results used in the present study. Transitic 
onset was determined by naphthalene flow visualizatioi 
The experiments of Muller and Bippes 23 and Deyhle an 
Bippes 22 employed swept flat plates under the action c 
favorable pressure gradients. The pressure coefficient mec 
sured in the experiment is shown in Figure 9 for the thre 
tunnels used in the tests. Transition onset was determine 
by the location where T = 0.5 

To facilitate comparisons with experiment, transitio 
onset was selected to be the point where T = 0.5. Thi 
assumes that the data reported in Radeztsky et al. 17 com 
sponds to a location where T = 0.5. As noted from Fig. 1C 
which compares present theory with the data of Radezt 
sky et al., 17 excellent agreement is noted for all Reynold 
numbers and surface finishes. For rms measured levels c 
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Table 3 Comparison of Re Vc , t r values at 50% intermittency predicted by the present method and measured in the 
swept plate experiments of Deyhle and Bippes. 22 


Facility 

Tu 

Plate Surface 

Experimental 
Re x .tr (xlO 5 ) 

Present 
Re x ,tr (xlO 5 ) 

Error (. %) 

NWB 

0.08 

Wooden plate. R z =■ 6^m 

6.5 

8.04 

23.7 

1MK 

0.15 

Plate with sandpaper. R z = 40 iim 

6.8 

6.12 

10.0 



Wooden plate,i2 2 = 6fim 

7.5 

7.49 

0.13 



Aluminum plate, sanded. R z = 5^m 

7.7 

7.63 

0.9 



Aluminum plate, polished. R z = 1.8 /im 

8.3 

8.48 

2.17 

1MK /screen 

0.27 

Aluminum plate, polished. R z = 1.8 iim 

7.8 

7.11 

8.85 

NWG 

0.57 

Aluminum plate, sanded. R z — 5am 

5.4 

2.18 

53.6 



Figure 9 Experimental pressure coefficient along the 
swept plate of Deyhle and Bippes 22 


Prwatnt (Q-2Stim AMS) 



Figure 10 Comparison of the present method with the 
experimental data of Radestsky et aL 17 


surface roughness, a sinusoidal distribution is assumed and 
a “peak-to-peak” roughness level is established by multi- 
plying the rms value by \/2. The calculatio ns assum e a 
freestream intensity, Tu = 0.09, reported by Dagenhart. 25 

Table 3 compares present predictions with measure- 
ments reported by Deyhle and Bippes 22 using three dif- 
ferent facilities and a variety of surface finishes. The 1MK 
tunnel is the 1 x .07m 2 DLR facility in Gottingen, the 
NWB tunnel is the 3.25 x 2.8 m 2 DLR facility in Braun- 


schweig, and the NWG tunnel is the 3 x 3 m 2 DLR facility 
in Gottingen. As may be seen from the table, excel- 
lent agreement is indicated for measurements taken in 
the 1MK facility. However, onset is overpredicted in the 
NWB tunnel and underpredicted in the NWG tunnel. The 
cause of the discrepancy in the NWB tunnel is not clear. 
On the other hand, the intensity in the NWG facility is 
high enough so that a by-pass mechanism may have been 
present. 

CONCLUDING REMARKS 

The present approach has developed, in the context of 
a CFD tool that employs turbulence modeling, a unified 
description for laminar, transitional, and turbulent flows. 
It allows for the presence of laminar fluctuations and treats 
them in a manner similar to that of turbulent flows. As 
a result, one can calculate the complete flowfield without 
having to use stability codes and at a cost comparable 
to that of existing CFD codes that employ two-equation 
turbulence models. 

Although the e n method or linear stability codes were 
not used in the calculations, results of stability theory 
played an important role in deter mining expressions for 
the eddy viscosity resulting from laminar fluctuations. Be- 
cause the physics underlying T-S or CF instabilities are 
different, corresponding stress-strain laws governing the 
flowfield are different. Further work is needed to develop a 
stress-strain law that encompasses the effects of all relevant 
instabilities. 
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